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Abstract. Let (L, h) — > (X, u) denote a polarized toric Kahler manifold. Fix a toric subman- 
ifold Y and denote by p t k '■ X — > R the partial density function corresponding to the partial 
Bergman kernel projecting smooth sections of L onto holomorphic sections of L k that vanish 
to order at least tk along Y, for fixed t > such that tk £ N. We prove the existence of a 
distributional expansion of ptk as k — > oo, including the identification of the coefficient of k n ~ 
as a distribution on X. This expansion is used to give a direct proof that if to has constant 
scalar curvature, then {X,L) must be slope semi-stable with respect to Y (cf. |RT 06 ). Sim- 
ilar results are also obtained for more general partial density functions. These results have 
analogous applications to the study of toric K-stability of toric varieties. 
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1. Introduction 

1.1. Polarized varieties and density functions. Let (X, L) be a smooth polarized projective 
variety of complex dimension n. Then the space = H°(X, &(L k )) of holomorphic sections 
of L k is a finite-dimensional vector space whose dimension grows like k n as k — > +oo. The 
ampleness of L also corresponds to the existence of a metric h on L of positive curvature, 
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Fh = —iu), where a; is a Kahler form. Denote by g the associated Riemannian metric. The 
choice of such a metric h with positive curvature equips Vk with a positive-definite innner 
product, making it a finite-dimensional Hilbert space. 

The density function pk : M — > R associated with (L, h) is defined by 

Pk(p) = ^2\e a , k (p)\ 2 , for pGM, (1.1) 

a 

where the {e a form an orthonormal basis of Vk- Here |e ai fc| 2 is the point-wise length-squared 
of the the section e a> k, computed using the metric h k on L k . We shall refer to |e a jt| 2 as the 
mass-density of e Qt k- It is easy to see that pk is independent of the choice of orthonormal basis 
{e a ,k}- 

The density function pk has been studied by many authors and it is known that for large k, 
there is a complete asymptotic expansion 

- £a,*-i, (1.2) 

J 3=0 

where the a,j are smooth functions on X which are local invariants of g, so for example 

do = 1 and a\ = -s, (1-3) 

where s is the scalar curvature of the Kahler metric g. The reader is referred to the literature for 
the background to these statements. Tian's famous paper [Tia90j essentially gives the /c n -term 



in this expansion; Zelditch [Zel98] obtained the complete asymptotic expansion (1.2) (see also 
|Cat99j . |MM07| and [BBS08J). The identification of a\ with half the scalar curvature appears 
in [LuOO] . 



The formula ( 1.2 ) is to be viewed as a local version of the Hirzebruch-Riemann-Roch formula: 
indeed, integration over X gives 

dim V k ~ A k n + A x k n - X + ■■■ , (1.4) 

where 



So 



A j= / a,—. (1.5) 
Jx n - 

A„ = (i)\o W = (1.6) 

Note that (1.4) also gives the leading coefficients of the Hilbert polynomial x{X,L k ) because 
the higher cohomology groups vanish for large k. The density function has been an important 
tool in the study of Kahler metrics of constant scalar curvature over the last decade or so, 
starting with Donaldson's pioneering work [DonOlj comparing balanced metrics — which make 
Pk constant for large k — with metrics of constant scalar curvature. The reader is referred also 
to |Don05b} IDon09| Biq06| IPS04[ IFinlO| for other contributions and aspects of this circle of 
ideas. 

In this paper we shall study a variant of pk , the partial density function associated in the first 
instance to a complex submanifold Y of X. Given a rational number t ^ 0, we may twist L k 
by JPy : , where J*y is the sheaf of functions vanishing along Y. This leads to a subspace 

V tk = H°(X, @{L k ) (8) J Y k ) (1.7) 

of Vk which of course inherits a Hilbert-space structure from Vk- Informally, Vtk is the space 
of sections of L k which vanish to order at least tk along Y. We shall obtain a distributional 
asymptotic expansion of ptk in the case that all data are toric and we shall use this expansion to 
give formula of the slope of the submanifold Y in the sense of Ross and Thomas [RT07| IRT06| — 
see below for definitions. From this we obtain an immediate proof that if the metric g in the 
Kahler class c\{L) can be chosen to have constant scalar curvature, then (X,L) must be slope 
semi-stable with respect to Y. 
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In order to state these results precisely, suppose that X is a smooth toric variety with cor- 
responding momentum polytope P C 1" and that Y is an irreducible smooth toric subvariety 
corresponding to a face F of P. The reader is referred to ^2] for a review of the correspondence 
between toric varieties and convex polytopes. We may choose affine coordinates (xj) on M. n so 
that 

F = {xi = ■ ■ ■ = x q = 0} n P, (1.8) 
while Xj ^ on P for j = 1, . . . , q. Let 

= XI H h x q . (1.9) 

Then <£> is non- negative on P and vanishes precisely on F. Pulled back to X, $ is everywhere 
non-negative, vanishing only on Y. In fact, $ vanishes quadratically on Y: if local complex 
coordinates z are chosen so that Y corresponds to z\ = ■ ■ ■ = z q = 0, $ = 0(Xa \ z i\ 2 ) ( see 



Lemma 2.4). 

Define three subsets of X: 

U t = ^[0,1), S t = $ _1 (i),A = S-^oo), (1.10) 

regarding <I> clS EL function on X. 

Denote by dat the Leray form of St, i.e. dat is an (n — l)-form satisfying datd^ = uj n /n\ 
along St- Let g be a toric Kahler metric on X with scalar curvature s. Then we may define a 
distribution a t on X (with support on St) by 

(a t J) = f^fdat- l -j t f s fm 2 g d<j t (f G C-(X)). (1.11) 

Then our first main theorem is as follows: 
Theorem 1.1. Zei f/ie notation be as above. Then 

Ptkip) = 0(k-°°) ifpeUt (1.12) 

and 

/Mp) = Pk(p) + 0(fc"°°) »/p g A (1-13) 
Moreover the O's are uniform if p moves in a compact subset respectively ofUt or Dt- 
If f £C°°(X), then 

- {L/ + k (I D , " + ( a -«) + h^- n) ■ < L14 > 

Rk is a torus-invariant distribution on X such that (R^,f) ^ C||/||c n+4 (x)! C being bounded 
independent ofk. 

Remark 1. In fact, our methods give a complete distributional asymptotic expansion of ptk, 
see Theorem 14.51 



1.2. Application to slope stability. The notion of slope stability of a polarized variety (X, L) 
with respect to a closed subscheme Z is due to Ross and Thomas [RT06. |RT07| . It was intro- 
duced as part of their study of K-stability of polarized varieties. Its advantage is that it may 
be relatively easy to show that (X, L) is slope-unstable with respect to a particular subscheme 
(or subvariety) implying that (X, L) cannot be K-stable. K-stability is important since it is 
conjectured to be the algebraic-geometric condition which is necessary and sufficient for the 
existence of a Kahler metric of constant scalar curvature in the Kahler class c\ (L) - at least if 
Aut(X,L) is discrete |Tia941 [Tia97l IDon02j . The necessity is known |Don05al [StoTffl] . but the 
sufficiency is a major open problem. 

One of our motivations for the study of the partial density function ptk was its application in 
the study of slope stability proposed by Richard Thomas and his coworkers |FKP + 09j . Theo- 
rem 1.3 realizes this proposal by giving a formula for the slope of a toric subvariety Y in terms 
of geometric data defined by the choice of a toric Kahler metric on X. From this formula, it is 
obvious that if g can be chosen to have constant scalar curvature, then X is slope (semi-)stable 
with respect to Y. 
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To describe these results more precisely, let us begin by recalling the relevant definitions. 
First of all, the slope fi(X) = fJ,(X, L) of a polarized smooth projective variety is defined in 
terms of the coefficients of the Hilbert polynomial 

X (X, L k ) = dim H l (X, G(L k )) = dim V k = A k n + A^ 1 + 0{k n - 2 ) (1.15) 

i 

as the quotient 

f i(X) = v(X,L) = ^-. (1.16) 



Note that by (1.6), we have 

v(X,L) = ±Av(s), (1.17) 

where Av(s) denotes the average scalar curvature on X. Alternatively, we can define n(X, L) 
by the formula 

dimH°(X,0(L k )) = A k n (l + n(X,L)k~ 1 + 0(AT 2 )) for k > 0. (1.18) 

Now let Z be a closed subscheme of X with ideal sheaf J?z- Then, if t ^ is such that 
tk G Z, we may consider the holomorphic Euler characteristic x(X,L k (g) )• This will be a 
polynomial of total degree n in the two variables k and tk and can therefore be written in the 
form 

X (X, L k ® ^f) = A (t)A; n + Axit)^' 1 + 0(^" 2 ), (1.19) 

where Ai{t) is a polynomial of degree at most n — i (and so is defined for all t £ 1). As 
previously, if t is fixed and k ^> 0, the higher cohomology groups vanish. So we also have 

dimV tk = A (t)k n + Ax(t)k n - 1 + 0(A^ 2 ) for fixed t, and A; > 0, (1.20) 

where we have written 

V tk = H°{X,0{L k )®jr z k ). (1.21) 
Since A)(£) and Ai(t) are defined for all real t, we may consider the quantity 

Mc(^) = Mc(^,i) = J ° f T;7., H 2 . (1-22) 

Jo ^o(i) dt 

provided the denominator is non-zero: this is called the slope of Z with respect to c. 

We can guarantee that Ao(t) > for t G [0,e(Z)), where e(Z) denotes the Seshadri constant 
of Z. Recall that one of the equivalent definitions of this quantity is 

e(Z) = sup {t : ir*L 0(-tE) is ample} , (1.23) 

where it : X — > X is the blow-up of X along Z and = ir~ 1 (Z) the exceptional divisor. 
The slope inequality for Z with respect to c is: 

fi c {jr z ,L)^fi(X,L). (1.24) 

Following [RT06], one makes the following 

Definition 1.2. (i) (X,L) is said to be slope semi-stable with respect to Z if the slope 

inequality ( |1.24[ ) holds for all c G (0,e(Z)]. 

(ii) (X, L) is said to be slope stable with respect to Z if we have strict inequality in (1.24) 
for all c G (0,e(Z)), and for c = e(Z) if e(Z) is rational and iI Q (X, G{L k ) ® ) 
saturates ^^ Z ^ k for k ^> 0. 

It is shown in [RTQ7] that, if (X, L) is K-semistable, then (X, L) is slope-semistable with 
respect to every closed subscheme Z and that, if (X, L) is (analytically) K-stable, then (X, L) 
is slope-stable with respect to every closed subscheme Z. In the light of the above-mentioned 
results (cscK implies K-stable, at least if Aut(A) is discrete) it follows that if there is a cscK 
metric in ci(L), then every closed subscheme Z of X is slope stable. 

Theorem 1 1 . 1 1 can be used to give a formula for /i c . Note further that, with this result, it is 
now easy to follow the approach of FKP + 09] to compute the difference /x c (JV> L) — fi(X, L): 
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Theorem 1.3. Let the data be as in Theorem \l.l\ Then we have, for < c < s(Y): 

f M c (jr Y ,L)- f i(X,L)= (2 J\ol(P t )dt^j j^ 6 (^( S -Av(s))d/i^ dt - * J \d<S>\ 2 g da c 

(1.25) 

In particular, if the scalar curvature is constant (s = Av(s)), we have 

H C (J? Y , L) - fx(X, L) < 0, for c G (0, e(Y)), (1.26) 
and (X,L) is slope semi-stable with respect to Y. 

Proof. The function 1 on X is a test function and so we may pair the distributional asymptotic 
expansion with it. The integrals of the local terms then give the coefficients Ao(t) and Ai(t): 

Mt)=(^) Vol(U t ), A 1 (t) = U^-\ (I adfig + AtY (1.27) 



where we have written At for (at, 1). Writing s = (s — Av(s)) + Av(s) and recalling (1.16), we 
obtain 



Now 



M*) = (J^f (voKUMX, L) + i ^ (s - Av(s))d/x 9 + 1 -A)j . (1.28) 

A t = I d^- 1 — / \d^\ 2 n da t (1.29) 
J St 2dtJ Si l '* { ' 



and by definition of dat, 



For e > 0, 

d. 
d7 



S t 



[ [ datdt = Vol(A) - Vol(i7 c ). 
Jo J s t 

\d$\ 2 do- s )jdt = \d$\ 2 g d(T c - ^ \d<S>\ 2 g da St . 



Now Id'&pdut is smooth and tends to zero near Y — its length with respect to g is \d<&\ g and <3? 
is quadratic in the distance toY. It follows that the integral over S e tends to zero with e since 
So = Y. The result follows by substitution of these formulae into (1.22). Indeed, the numerator 



of ( |1.22[ ) is 

M »(X,L) j\o\{U t )dt 



A 1 (t)dt+±(A (c)-A (0))= 2 _ 



1 (\_ 
+ 2 V 2vr 



(1.30) 



Dividing by f c Ao(t)dt as given in (1.27) and rearranging, we obtain the formula stated in the 
the theorem. We obtain the slope semi-stability of (X, L) with respect to Y from the strict 
inequality ( |1.26[ ) by noting that c 1— > /u c (j^y , L) is continuous for c € (0, e(^)]- □ 

1.3. More general partial density functions. There is a generalization of Theorem |1.1| to 
partial density functions associated with more general subspaces Vtk of We shall describe 
the set-up we have in mind in the toric case. 

Let X correspond to the polytope P as above, and suppose that we have 1-parameter family 
of subpolytopes Pt of P. More precisely, suppose that 

P t = Pn p|{$ a (s) >t) (1.31) 

where the <3? a are affine-linear functions on K n with rational coefficients and A is some finite 
index set. 

We assume that Pq = P and are interested in values of t for which Pt is an n-dimensional 
polytope strictly smaller than P. 
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Because the coefficients of the <& a are rational, there is a positive integer N such that kPt is 
an integral polytope for all positive integers k divisible by N. For such k, it therefore makes 
sense to consider the subspace Vtk of holomorphic sections of L k corresponding to the points of 



[k 1 'Z) n r\Pt. This subspace defines a partial density function ptk as before. Theorem 5.2 below 



extends Theorem [TTT] to this more general setting. We defer the precise statement, but note here 



that the picture is similar to the one we see in Theorem 1.1 there is a tubular neighbourhood 
D t of a reducible subvariety of X, with boundary St, such that ptk(x) is rapidly decreasing in k 
if x G Dt, ptk(x) is essentially equal to pk(x) if x Dt, and there is a distributional expansion 
of ptk with an explicit contribution at supported on St- The set St is now no longer smooth, 
however (it has singularities at points mapping to intersection points of two or more of the 
hyperplanes $ a (x) = t) and there is an additional term in at which is supported on the singular 
locus of St- 



Just as Theorem 1.1 gives information about slope-stability if the metric is of constant scalar 



curvature, so Theorem 5.2 gives information about K-stability. Indeed, an analogue of Theo- 
rem 1.3 is Theorem 5.6 which gives a formula for the Donaldson-Futaki invariant of a toric test 
configuration in terms of the distributional expansion of ptk ■ This formula immediately implies 
that if the metric has constant scalar curvature, then {X, L) is K-polystable with respect to 
every toric test configuration. This result was previously proved in [ZZ08J without the use of 
density functions. 



Remark 2. Theorem 5.2 can also be used to give a formula for the slope p c (j?z,L) for more 
general ideals J*z- We have omitted an explicit treatment of this, however, because the appli- 
cation to K-stability seems to be more interesting. 

1.4. Relation to previous work. Asymptotic expansions of what we are calling partial den- 
sity functions were studied in detail by Shiffman and Zelditch [SZ04]. Their point of view was 
that of random polynomials with prescribed Newton polytope, and the partial density functions 
then appear as 'conditional expectations'. Our results on the distributional expansion of ptk 
go beyond those of |SZ04| . by giving information about p t k at the interface St between Ut and 
Dt- On the other hand, the results of Shiffman and Zelditch in the interior of these regions 
are much more precise than ours. We mention also that these authors deal only with the case 
that X = CP" with the Fubini-Study metric (though the extension to general toric metrics is 
probably straightforward) and that their methods are completely different from ours, the start- 
ing point being the description of the Szego kernel as a Fourier integral operator with complex 
phase. By contrast, our methods are elementary and explicit. 

Moving away from the toric case, Berman |Ber09| announced that in general, given a complex 
submanifold Y C X and with Vtk defined as in ( |1.7[ ), there exist open subsets D t and Ut of X 
satisfying the conditions (4.16) and (1.12) of Theorem |1.1[ However, in this generality, there is 



no information about the smoothness of dDt nor about the 'transition behaviour' of ptk near 
dD t - 

Finally we note that this work grew out of the first author's Edinburgh PhD thesis [Pokll] 
which contains further pointwise information about the asymptotic expansion of toric partial 
density functions. 

1.5. Outline. The remainder of this paper is organised as follows. In ^2j we collect some 



standard notions from toric geometry. The key to our subsequent analysis is the formula (2.14), 
which expresses the mass-density of a unit-norm basis element e a ± of L k in terms of a function 
<p(at, y) derived in simple fashion from the symplectic potential u which defines our Kahler 
structure. 



In we use Laplace's method to compute a distributional asymptotic expansion of \e a 



2 

fe > 



following closely the approach of [BGU10J. The key results here are Propositions 3.3, 3.4, and 
3.5 In Q these results are combined with the Euler-Maclaurin formula for (lattice) polytopes 



to obtain Theorem ll.il 

In ^5] we shift attention to the more general partial density functions mentioned in §L3j The 
method used to obtain the distributional asymptotic expansion in this case is the same as that 
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followed in £[4j it is, however, technically more complicated to obtain a nice formula for the 
distributional term at in this more general case. The remainder of the paper is devoted to the 
application of Theorem 5.2 to obtain a formula for the Donaldson-Futaki invariant of a toric 
test configuration and to deduce that cscK implies K-polystable with respect to such toric test 
configurations. 



1.6. Acknowledgement. We thank Julius Ross, Richard Thomas and Steve Zelditch for useful 
conversations. The second author was supported by a Leverhulme Research Fellowship while 
this work was being completed. 



2. Background 

We review very briefly the elements we shall need of toric geometry, referring the reader to 
|Ful931 1(;ui941 IAbr981 IAbr03| IBGUlOj for more details. 

2.1. Combinatorial description of toric varieties. First of all, we recall the correspondence 
between smooth projective polarized toric varieties (X, L) on the one hand and integral Delzant 
polytopes on the other. Thus we suppose that X is a smooth (connected) projective variety of 
complex dimension n and that X contains a dense open subset X° isomorphic to a complex n- 
torus T™ . We suppose further that the standard action of T™ on itself extends to a holomorphic 
action of T™ on I. 

Now suppose that T n C T™ is a compact real torus. We are interested in Kahler structures 
on X that are invariant under T n , so that T n acts by isometries of the Kahler metric g. As the 
action is holomorphic, the Kahler form cj is then automatically T n -invariant. 

In this setting there is a moment map fx : X — )■ t* (where t* is the dual of the Lie algebra 
of T n which is isomorphic to M n ), and the image P = fi(X) is a convex compact polytope, the 
convex hull of the images of the fixed-points of the T n -action. The restriction of /x to X° is a 
fibration with image the interior Int(P) of P, with fibre T n . 

The Lie algebra t of T contains the weight lattice A = Ker(exp) and correspondingly t* 
contains the coweight lattice A*. 

The condition that X is smooth and has a given polarization L translates into the condition 
that P is an integral Delzant polytope: 

Definition 2.1. A convex polytope P C t* is Delzant if 

(1) There are n edges meeting in each vertex v. 

(2) The edges meeting in the vertex v are rational; i.e., each edge is of the form v + tei, with 
t^ 0, t G E and a G A*. 

(3) The e%, . . . , e n in (2) can be chosen to form a basis of A*. 

An integral Delzant polytope in t* is a Delzant polytope whose vertices lie in A*. 

If P is a Delzant polytope, we may write P as the intersection H\ n • • -dH^ of a finite number 
of affine half-spaces and we may assume that for each a, 

H a = {x G t* : £ a (x) := (x, v a ) - \ a > 0}, 

where v a G A is primitive. The intersection of the boundary of H a with P defines a codimension- 
1 face or facet of P, which will be denoted Q a : 

Q a = {x £ i* : £ a (x) = 0} H P. (2.1) 

P is integral if and only if all the A a are integers. 

More generally, if Q is any face of codimension q of P, there will be a subset {oi, . . . , a q } C 
{1, . . . , a 1 } such that 

Q = Q ai n---nQ aq . (2.2) 

Note that the conormal space A^Q (that is, the annihilator in T*M n of TQ) is just the span of 
{v ai ,...,u aq } (or equivalently of the {d£ ai d£ aq }. 
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2.1.1. Leray forms. 

Definition 2.2. Let Q a be a facet of P. The Leray form da a of Q a is the (n — l)-form on 
Q a with the property that da a d£ a = dx (Lebesgue measure) on Q a . Let Q a b = Q a n Qb be a 
co dimension- 2 face of P. The Leray form dr a & of Q a b is similarly defined to be the (n — 2)-form 
on Q a b such that dr^d^d^ = da;. 

In order to keep our notation short, we shall denote by da the measure on dP whose restriction 
to the relative interior of Q a is da a and by v the almost-every where defined section of T*R n 
such that v is the conormal to Q a on its relative interior. The measure dr with support on the 
(n — 2)-skeleton of dP is defined in the analogous way. 

In we shall need to consider polytopes which are not simple (so that more than n facets 
can come together in a vertex). Note that for any convex polytope, however, every codimension- 
2 face is always the intersection of just 2 facets and so the Leray form dr is still well-defined 
in this case. If the polytope is simple, then every face of codimension q is the intersection of 
precisely q facets, and so has a well-defined Leray form. 

2.1.2. Adapted coordinates. 

Definition 2.3. If p G P, there is a unique face F of P which contains p in its relative interior. 
This relative interior will be denoted by F p . 

In particular, p £ F p . The two extreme cases are F p = Int(P) if p is an interior point of P, 
and F p = {p} if p is a vertex of P. 

If p E P, adapted coordinates centred at p will mean a choice of affine coordinates x on W 1 
such that 

• F p is an open subset of {x\ = X2 = ■ ■ ■ = x q = 0}; 

• Xj ^ on P if j = 1, . . . , q; 

• the point p corresponds to x = 0; 

• Lebesgue measure on W 1 is given by dx\ • • • dx n . 



It is clear that such coordinates always exist: if F p is the relative interior of the face Q in (2.2), 
then we take Xj = t a . for j = 1 , . . . q and choose the remaining coordinates so as to satisfy the 
remaining conditions. Then the Leray form of F p is just dx g +i • • • da; n . In the case of a Delzant 
polytope, these coordinates can be chosen so that A* is identified with the standard lattice Z n , 
(in other words so that the change of coordinates lies in SL(n,Z)). The fact that P can be 
covered by a finite system of adapted coordinate charts will be useful in the next section. 

Note that an adapted coordinate chart gives rise to a smooth system of (local) coordinates 
on X in the following way. Letting (8i, . . . ,9 n ) be angular variables dual to the coordinates 
(xi, . . . , x n ) (i.e. the 9j give coordinates on t), then the real and imaginary parts of ^Jx~je %e i , for 
j = 1, . . . , q extend to be smooth near /i" 1 (F p ); and (xj,6j) for j = q + 1, . . . , n also lift to be 
smooth functions near /x _1 (F p ). 

Hence we have the following result: 

Lemma 2.4. Let P be an integral Delzant polytope as before and let £ be the defining function 

of a facet Q. Then fJ,*(£) is smooth on X and vanishes quadratically on Y = ^ 1 {Q) and is 

positive elsewhere on X. 

More generally, if a face Q of P is defined by the subspace x\ = ■ ■ ■ = x q = 0, the Xj being 
elsewhere on P, then n*(x\ + • • • + x q ) is smooth on X , vanishes quadratically on fi~ 1 (Q) 

and is positive elsewhere on X . 

2.1.3. Lattice points and holomorphic sections. Two Delzant polytopes P and P' determine 
isomorphic toric varieties if they are combinatorially the same and the set of normals to the 
facets of P is the same as the set of normals to the facets of P'. The polytope itself fixes in 
addition a Kahler cohomology class [up] £ H 2 (X, M) which is integral if and only if the polytope 
is integraQ In this case, there is a T^-invariant holomorphic line bundle L = Lp on X such 



"*"To be more precise, we should say that [cup] is integral iff fi (which is only determined up to an additive constant) 
can be chosen to make P integral. 
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that ci(L) = [ojp] and it is well known that the T"-equivariant sections of L are in one-one 
correspondence with the points of PPi A* and form a basis of H°(X, 0(L)). Replacing L by L k 
corresponds to replacing A* by the rescaled lattice 

A* k = {y G t* : ky G A*} (2.3) 

so that there is a basis of sections s a ^ of H°(X, &{L k )) indexed by a G PC\K* k . (Alternatively, 
we can think of the lattice as fixed and replace the polytope P by the dilated polytope kP to 
get this basis of sections.) 
It is worth recalling that 

s a ,k(y) ^O^aePj (2.4) 

and that the morphism X — > FH°(X, &(L k )) defined by this basis of holomorphic sections is 
an embedding. 

Let F be a face of P, of codimension q. Then Y = /z _1 (.F) is a toric subvariety of X, and 
conversely any irreducible toric subvariety of X is equal to n~ l (F) for some face Q. If F is 



written as in (2.2), then the normals u ai , . . . , v a G A generate a subtorus Tp of T n . This is the 



stabilizer of Y which is toric with respect to the quotient torus T n /Tq. 

2.1 A. Seshadri constant. If Y is the subvariety corresponding to the face F of P, with F 
defined as usual by the condition <3?(x) := x\ + ■ ■ ■ + x q = in adapted coordinates, then from 
the definition ( 1.23[ ), the Seshadri constant e(Y) of Y is given by 

e(Y) = sup{t > : <&(x) = t contains no vertex of P}. (2-5) 

2.2. Toric Kahler metrics. A choice of toric Kahler structure on X corresponds to choosing 
a symplectic potential u on P (see |Abr 03j ) . Thus u : P — > M is a strictly convex function, 
smooth in the interior and satisfying the boundary condition 

1 d 

u(x)--J2U^og£ a (x)eC°°(P) (2.6) 

a=l 

(i.e. this difference is smooth up to the boundary of P). Given such a symplectic potential, set 

H = Hess(u), in other words Hij = didjU, (2-7) 

and 

G = H~ 1 . (2.8) 
In terms of these matrices, the Kahler structure over X° = ^ _1 (Int(P)) is given by 

u) = dxj A d6j, g = Hijdxidxj + G ij d9id6j, (2.9) 



and the boundary condition (2.6) ensures that this extends smoothly to X. 



Although u is not smooth up to the boundary of P, the restriction up of u to any face 



F of P is well-defined by (2.6). As part of the condition of 'strict convexity', up is required 
to be strictly convex and smooth in the interior of F and to satisfy the analogous boundary 
conditions. In fact, up is the symplectic potential for the restriction of the Kahler structure to 
the toric submanifold n~ l (F) of X. 
The function 

1 d 

u (x) = -^4Wlog4(x) (2.10) 

0=1 

is strictly convex in P and clearly satisfies ( |2.6[ ) . This symplectic potential gives a special choice 
of toric Kahler structure on X called the Guillemin metric on X [Gui94j. 

Note that the addition of an affme-linear function of x to u does not affect the metric. It does, 



however, affect the metric on the line bundle whose curvature is the Kahler structure (2.9) 



Definition 2.5. Denote by s a ^ a choice of section corresponding to the lattice point a, nor- 
malized so that the maximum value of \s a ^(y)\ 2 is equal to 1. 
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For each a and k, s a: k is thus defined up to multiplication by a unit complex number. 
Following [BGUlOj . define 

ip:Pxlnt(P) — > R, tp(x,y) = 2(u(x)-u(y)-(du{y),x-y)). (2.11) 

Then we have the key formula [BGU101 ISZ07] 

\saM\ 2 = e~ kv{a ' y) (2-12) 



for a section s a k normalized according to Definition 2.5 

We note that for fixed y, x i— > tp(x, y) differs from u{x) by an affine function of x. In particular, 
it is strictly convex. We also have 

<p(x, x) = 0, V x (p(x, y) = if x = y and V y tp(x, y) = if y = x. (2-13) 

It follows from the convexity in x that tp(x,y) with equality if and only if x = y, at least 
for y E Int(P). 

The s a) k are automatically mutually orthogonal with respect to the 1? inner product, and so 
rescaling by the length of s a ^ we obtain an L 2 -orthonormal basis of sections e a &, satisfying 

e -k(p(a,y) 

l6a ' k{y)l = (2^/ p e-M^)dz- (2,14) 



We shall refer to \e a k\ 2 as the mass-density of e, 



Remark 3. The formula (2.12) continues to be valid, with a suitable extension of the definition 



of tp, when y E dP. This requires some care: indeed, we see that for (2.4) and (2.12) to remain 
consistent, we need to define tp(x, y) = +oo if a F y . 

3. Asymptotic expansion of the mass-density function 
The goal of this section is to obtain the large-/c asymptotic expansion of the quantity 

|2 Jpe-M«,»?/(y)dy 

{l6aM > J) - f pe -k V (a,y) dy ' 



where e a ^ and tp are as in (3.1) and (2.14) and / is any smooth function on P. 

We shall use Laplace's method for this, but this entails an understanding of the critical 
points and some other global properties of the function y \— )■ tp(a, y) for fixed a. The analysis 
is straightforward if a is an interior point of P but a bit more complicated if a lies on the 
boundary. 

We follow the argument of |BGU10| closely here. We have nonetheless provided the details, 
because their discussion applies only to the Guillemin metric on X; and on the other hand 
Sena-Dias fSDIOj provided the extension to general toric metrics but did not fully analyze the 
situation at the boundary. 

The following will be used in this section (and the rest of the paper): 

• H is the Hessian didjU of u, and G = H ; 

• a Euclidean structure is fixed on W 1 , the length of a vector v being denoted by \v\; 

• we denote by ||/|| r the C r -norm defined by our given Euclidean structure. 

3.1. Properties of u and tp. We begin with a statement of the properties of H and G that 
will be needed later. 

Lemma 3.1. (i) There is a constant c > such that 

(H(x)v, v) ^ c\v\ 2 for allx£P,v£ R n , (3.2) 

where the LHS has to be interpreted as +oo if v TF X ; 
(ii) G = H~ 1 is smooth on P, G{x) is positive-semidefinite for all x E P and 

G(x)£ = if and only if £ E N*F X (3.3) 

(i.e. £ is conormal to F x at x.) 
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Proof. Let p G P and choose adapted coordinates such that F p is defined by the vanishing of 
xi,...,Xg. In particular these functions are ^ on P. It is convenient to write x' = (x\, . . . , x q ) 
and x" = (x q+ i, ...,x n ). 
Set 

D = diag(xi, . . . ,x g ). (3.4) 
Then corresponding to the splitting of variables x = (x',x"), we have the block decomposition 

*-(««>- (<% + * £) (3-5) 

of the Hessian of u, where 

'h\ hI) < 3 - 6 > 



is smooth. At the boundary, H2 is the Hessian of up , H2 is positive-definite near F p (cf. {2.2). 
Hence H is positive-definite and (H(x)v, v) = +00 if and only if v has a non-zero component in 
the subspace spanned by ei, . . . , e q , i.e if v is noi tangent to F p . Covering P by a finite number 
of open sets of this kind, a simple compactness argument establishes part (i) of the lemma. 
For part (ii), let 

[V2D 1 /* \ 



A 
Then 



V H?l\ 



AHA = 1 + R(D) (3.7) 

where 

/ 2D^H D^ V2D^ Hl H-^\ 
R{D) -[V2H 2 ^HiD^ J' ^ 

Now certainly = 0(|x'| 1//2 ) for small x' and so sufficiently close to F p , we have 

00 

(1 + R(D))- 1 = Y,(-R(D)) J = 1 + S(D), (3.9) 

j=0 

say. It is easy to see, moreover, that 

A- 1 S(D)A- 1 = f)s(D)(° ^ (3.10) 
where S(D) is now a smooth function of x' . It follows that the inverse G of H has the form 

2 „%"0 + (o ?H(o J). (3.U) 

In particular G is smooth up to the boundary and everywhere positive-semidefinite. 

For the last part, suppose first that Q is a facet of P and suppose also that coordinates are 
chosen so that Q = {x\ = 0}. Now let p £ Q. Then F p C Q and so x\ will be among the 
coordinates adapted to F p and centred at p. With these choices, if £ annihilates TQ then it 



must be a multiple of e±, and by (3.11), 



Ge 1 = 2x l e 1 + x x r^ fj S(D)e 1: (3.12) 

which shows that Ge\ = at p. It follows that Ge\ = on the whole of Q (since p £ Q was 
arbitrary) . 

If now F = Q\ n • • • n Q q is an arbitrary face of P, then choosing adapted coordinates, we 
know that Gej = along Qj, and so Gej = for all j = 1, . . . q on F. Since N*F is the span of 



{ei, . . . ,e q }, the 'if part of (3.3) follows 



The 'only if part of (3.3) is proved similarly. □ 



We now give some key properties of (p. 
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Lemma 3.2. (i) The function ip(x,y) is smooth on P x Int(P) and there is a constant 
c > such that 

ip(x, y) ^ c\x — y\ 2 for all x G P, y G Int(P). (3.13) 

(ii) The function ip extends naturally to a function on P x P with values in [0, oo] such that 
(p(x,y) = +oo if and only if x G" F y and satisfying 

<p(x, y) ^ c\x — y\ 2 for all x,y G P. (3-14) 

(Hi) Let p be a point of P and let x = (x',x") be adapted coordinates centred at p. Then 
there is a constant C such that for all sufficiently small y, 

<? 

3=1 

in these coordinates. 

Proof. Suppose first that y G Int(P). Let v G W 1 be any unit (with respect to our arbitrary 
Euclidean structure) vector, and define 

f(t) = <p(y + tv,y). (3.16) 



The domain of / is the interval / such that y + tv G P. In particular, G / and by (2.13) 

/(0) = 0, f'(0) = 0, (3.17) 

and 

f"(t) = {H(y + tv)v,v)>c (3.18) 
by Lemma 3.1 Integrating this from to t and using Q3.17D , we obtain f(t) ^ ct 2 /2. Since 
\x — y\ = \t\, part (i) follows. 

To understand the behaviour of (p near the boundary, let us write 

u = uq + w, (3.19) 



where no is the Guillemin potential (2.10) and w is smooth on P. A simple computation gives 
<p(x,y) = J2(U^og£ a (x)-log£ a (y))-£ a (x-y)) + ^(x,y), (3.20) 

where 

ip(x, y) = 2[w(x) - w(y) - (Vw(y),x - y}]. (3.21) 
is smooth, hence bounded, on P x P. If we fix x and let y — > yo G dP, then it is clear that 
p(x,y) — > +oo if there is an index a with £ a (yo) = but £ a (x) > 0. This is precisely the 
condition x F yo which is consistent with (2.4) and (2.12). 

It remains only to consider the situation that x G F y where F y is the interior of a proper face 
of P. Now the restriction <pp, say, of (p to F, is given by 

p F {x,y) = 2(u F (x)-u F (y)-(Vu F (y),x-y)), [x G F,y G Int(F)). (3.22) 

where u F , the restriction of u to F, is the symplectic potential for the restriction of the Kahler 
structure to u _1 (F). 

Thus we can replace P by F, u by u F and <p by ip F in the argument at the beginning of this 
proof to obtain (3.14) for x G F y . This completes the proof of part (ii). 
The last part is a local computation. In adapted coordinates, 

q 

<p(0, y) = 2[«(0) - u(y) + (Vu(y),y)] = ]T y, - 2{v(y) - v(0) - (Vv(y), y)}. (3.23) 

i=i 

Now the part in curly brackets is a smooth function of y which vanishes and has gradient at 
y = 0. Hence for sufficiently small y, we can bound this by a multiple of \y\ 2 , giving 

<? 

^(0,yK]Ty,+C|y| 2 
i=i 

as required. □ 
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3.2. Distributional asymptotic expansion of |e aj fc| 2 . The main goal of this section is the 
following: 

Proposition 3.3. Let f be a smooth T n -invariant function on X and denote by the same letter 
the corresponding function on P. Denote by s the scalar curvature of the metric with symplectic 
potential u. For each a G P n A* k , recall that e a ^ is the unit-length holomorphic section of L k 
corresponding to the point a. Then we have 



K, k \ 2 J) = /(«) + ^ (*(«)/(«) + IWGPfXa)) + ^(Rk(a)J) 



(3.24) 



2k \~ x ~~ yJ v ~' ' 2~'~ JV ~ " /v "7 ' k 2 
where Rk(ct) is a distribution which satisfies 
. for fixed a G P, (R k (a)J) ^ C\\f\\c*; 

• for each fixed test-function f , (R k (a), f) is smooth in a and bounded for > 0. 

Here we recall that G = {G^) is the inverse of the Hessian of the symplectic potential and 
that s is the scalar curvature of the metric g. 

We begin with a stronger result covering the case that supp(/) does not contain a. 



Proposition 3.4. Suppose that a £ P and f G C°°(P) with a G" supp(/). Then 

(\e a , k \ 2 ,f)=0(k-°°) 

for large k. 
Proof. Since 



i\e a ,k\ 2 J) 



/ e -M«.w)/(j,)dj/ 



(3.25) 



(3.26) 



f p e-M<*,2/) dy ' 

we need an 0(k~°°) upper bound for the numerator and an 0(k N ) lower bound for the denom- 
inator. In fact we shall obtain an exponentially small upper bound for the denominator. 



By Lemma 3.2 



\e~ k ^ z) f{z)\ < e - ck \ a - z \ suppl/l 
and so if the distance from a to supp(/) is d, 



- kip{a ' y) f(y)dy 



-ck\x~y\ 2 i 



\x-y\^d 



Now 



/ ■ 

J\z\^d 



- fcc|z|2 dz € Ce~ kd2 



for some constant C independent of k. Hence 

e -M«.*)/( 2 ) dz 



< Ce 



-kcd 1 



(3.27) 
(3.28) 
(3.29) 

(3.30) 



We complete the proof by obtaining a suitable lower bound on the total mass of s a ^. Suppose 
that y = (y',y") are adapted coordinates centred at a so that we are in the situation of part 
(iii) of Lemma 3.2 Suppose further that the subset 

V = {0 < Vj < e for j = 1, . . . ,q} x {\y"\ < e} (3.31) 

(3.32) 



is contained in P. Now by (3.15), 



3=1 

and by shrinking e if necessary, we may absorb the |y'| 2 -term into the linear term, getting 

i 

<p(0,y)^2^y j + C\i/ , \ 2 (3.33) 

3=1 
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for y € V. Then 



jf e-fc^dy^ J v exp(-k\2^y j + C\y''\ 2 \\ dy. (3.34) 

Now the difference between this integral and the integral over x ~S. n ~ q is exponentially small 
in k, with a constant depending upon e). Since 



o 



e-^^and / e "*V = ( 3 .35) 



it follows that 

' e -k V (a,y) dy ^ Ck -q-(n-q)/2 = ^-(n+g)^ (3.36) 



Dividing (3.30) by (3.36) completes the proof. □ 

The effect of this Proposition is to localize this (\e a: k\ 2 ,f) to an integral over an arbitrar- 
ily small neighbourhood of a in P, up to exponentially small terms. We now calculate this 
contribution recursively. 

With a fixed as before, choose adapted coordinates y = (y',y") as in the previous proof, and 
let V be as in (3.31 ). Choose a cut-off function p G Cq°(V), ^ p ^ 1, where p = 1 in a smaller 
neighbourhood W = \V of a in P. 

Define the operator 5 : Cg°(F) -»■ C°°(y,IR n ) 

<WG/)= f djf(ty)dt (3.37) 



so that 

/(y) - /(0) = VjSjfiy) (3.38) 
(summation convention) for all y E V. For any function / 6 C^°(y) define the linear operator 
: Cg°(F) -> Cg°(V) by 

&f = \d i (pG i ^6 j f). (3.39) 

Note that the operator depends also on the point a. When we need to draw attention to this 
fact, we shall denote it also by 3i a . 

The significance of this operator is as follows 

Proposition 3.5. With the notation as above, we have, for any N ^ 1, 

N 

(\e a ,k\ 2 , /) = E fc " m ^C/(«) + k-^&N+itM, (3.40) 

m=0 

where the remainder term &N+l,k,a(f) is smooth in a for fixed f and satisfies 

^N+l,k,a(f)<C N \\f\\ C 2N+2. (3.41) 

uniformly in a and k. 

Proof. Given the test-function /, write 

f(y) = p(y)/(0) + P (y){f(y) - /(o)) + (l - P (y))f(y), (3.42) 

and substitute this into j e~ klfi( - a ' y ^ f(y) dy, getting 

e -^(y)f(y)dy = /(0) J e- k ^p(y) dy + J e~ k ^ p(y)(f(y) - /(0)) dy (3.43) 

+ Je- k *to(l-p(y))f(y)dy. (3.44) 



In the second term, use (3.38) and note also that 

Q. e -Wv) = -2kH ijyi e~ k ^ , so that ^G ij dje'^ = -%e"^ (y) . (3.45) 
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Hence 



" kv(y) p(y)(f(y) - /(0)) dy = -k- 1 J ±G*a i e- k rp6 j fdy = k' 1 J e^&fdy, (3.46) 



where we have neglected the boundary term 



G ij is je - k(fi P 5ifda. 



dp 



(3.47) 



This is justified because G l3 Vj = on the interior of each facet of P (see part (ii) of Lemma 3.1 ). 
In summary, then, we have the formula 

-^(a,y) f{y)dy = f{a) J e -^(a,y) p{y)dy + k -l J e -k<p(a,y) ^ f {y) dy 

+ I e- k ^y\l-p{y))dy. (3.48) 



We can now iterate: we apply (3.48) to the second term on the right-hand side, (i.e. with f(y) 

e -^(a,y)^N+l f{y)dy 



replaced by @ a f(y)). After ./V steps, we obtain the formula 

e- k ^f(y)dy = ^ N , a (f)(a) [ e~ k ^ p{y) dy + AT*" 1 



+ 



- k ^ a ' y \l-p(y))^ N ,af(y)dy. 



From the proof of Proposition |3.4 we have 

e- k ^y\l-p{y))dy = e- ck * m {a) 



(3.49) 



(3.50) 



for some c > where rjk(a) is smooth in a and uniformly bounded in k provided that a moves 
in some smaller subset \V , say. Moreover, (3.34) and (3.35) imply that 



(^j e -k<p(a,y) dy ^ J e -k<p(a,y) dy = 1 + 



(3.51) 



where rf k has the same properties as 

Hence, dividing by f e~ k{p , we get ( 3.40[ ), where 

&N+i,kAf) = k N+1 Vk(a)e- ck ^ N , a f(a) 

+ (J e -Ma,y)\ J e -Ma.») {^f(y) + k N+1 (l - p(y))^ a f(y)} dy. 

(3.52) 

It is clear from this formula that for fixed / and k, &N+i,k,a is a distribution in / depending 
smoothly on a. 



On the other hand, by (3.36) and (3.50) we have 

\&N+i,kAf)\ < C(sup|< +1 /| +e- cfc£ sup|^v, Q /|) 



(3.53) 



directly from (3.52). Since the operator 5 has the same boundedness properties as a differential 
operator, 

[|*/Hr < AWfWcr+x (3.54) 
for r ^ 0, where A = A r is some constant, it follows that the operator *2> behaves like a 
second-order operator in the sense that we have an estimate: 

Wllr < A\\j\\ C r + * (3.55) 

(for some different constant A = A r ). It follows by induction that sup \ S>™f\ is bounded by a 
multiple of ||/||c2m+2. The estimate 

\\^N+l,k, a f\ < C\\f\\ C 2M + 2 (3.56) 

now follows by combining these observations with ( |3.53 ). □ 
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To obtain Proposition 3.3 from this expansion, we take N = 1, getting 

{K,k\ 2 , f) = /(«) + k- x 9f{a) + k- 2 <% 2 (f), (3.57) 

and it follows from the formula for &2 that this error term has the stated properties. It remains 
to compute 3>f(a). In local coordinates, with a corresponding to as before, 

Sjf(y) = 3,-/(0) + ^ WW + 0(|y| 2 ) (3.58) 
from the Taylor expansion of f(y) — /(0) and, after a little manipulation, we obtain 

0/(0) = ^(/G«)(0) - 1/(0)^(^(0). (3.59) 



The formula (3.24) now follows by inserting Abreu's famous formula for the scalar curvature, 

1 

S = ~2 



-didjG tj . (3.60) 



4. Proof of Theorem 11.11 

We now bring the ideas of the previous sections together to prove Theorem |1.1| Recall that 
the setting for that Theorem was as follows: 

• A toric variety X with moment polytope P; 

• A face F = Q\C\ ■ ■ ■ Q q with Qj defined by Xj = for j = 1, . . . , q. 

• The subpolytope Pt = P D {<&(x) ^ t}, where $(x) = x% + • • • + x g . 

Then our partial density function is given by 

ptk(y)= Yl \ e ^(y)\ 2 

oeP t nA* 

(regarded, by abuse of notation, as a function of y G P), where the terms in the sum are given 
by gjg. 
Define 

C t = ^[O,*),^*) = $-Ht),P t = Q^faoo). (4.2) 



These are the subsets of P corresponding respectively to the three subsets Ut, St and D t in ( 1.10 ). 
By torus-invariance, it is clearly enough to prove the 'pushed-down' version of Theorem 1.1 i.e. 
to work entirely on P. 

We begin by establishing the first part of the Theorem |1.1[ namely the equations (1.12) and 
(1.13) restated as follows: 

Proposition 4.1. Let K be any compact subset of Ct- Then 

Ptk( x ) = 0(k~°°) uniformly for x G K (4.3) 
and if K' is a compact subset of Pt, then 

Ptk{x) = p(x) + 0(k~°°) uniformly for x € K (4.4) 
Proof. If x G K and a G Pt, we have 

\e a , k \ 2 ^Ce- kd ^\ (4.5) 

where d denotes Euclidean distance. Summing over lattice points of Dt gives the result, since 
the number of lattice points is 0{k n ). The proof of the other part is the same, the roles of Ct 
and Dt being interchanged. □ 

Remark 4. If X = CP n with the Fubini-Study metric, then more precise pointwise estimates 
of this kind are given in |SZ04j . at least for points x in the interior of P. There, D[ is called the 
'forbidden region'. 
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4.1. The Euler Maclaurin formula. In order to obtain an expansion in powers of k from 
(4.1 ), we use the Euler-Maclaurin formula to replace the sum over lattice points by an integral, 
up to a controlled error term. The version we use is as follows: 

Theorem 4.2. Let P be a convex integral polytope of dimension n, with integral conormals. 
Let Lebesgue measure dx be normalized so that the integral of the unit cube in Z n has volume 1 
and let da stand for the Leray form of dP. Then we have 

Y f(a) = k n / f(x)dx + — / fda + k n ~ 2 E k (PJ) (4.6) 
PnA k Jp 2 JdP 

where E k (P, f) is bounded by a multiple of Vol(P)||/||2n (the C 2n -norm of f again). 

Remark 5. Since many results of this kind are available in the recent literature, we shall 
be content to sketch a proof. Following the method used by Donaldson in the appendix of 
|Don02j . we reduce to the case that P is a lattice simplex. Then we are content to quote the 
Euler-Maclaurin formula with remainder from [KSW03J to complete the proof. 

We note references such as [GS071 IKSW03] ) give complete asymptotic expansions of lattice 
sums at least if P is a simple polytope. The theorem stated here applies to any lattice polytope, 
and this will be important in This simple statement (4.6) should be viewed as an extension 
of the 'trapezium rule' (with remainder) for approximate integration of functions of one variable. 

Proof. Given / E C°°(P), consider on the one hand 

S k (f,P) = £ f(a)- 1 - Y, /(«)■ 

aePnA fe ae9PnA fc 



I(f,P) = j p f(x)dx. 



and 

We aim to show first that 

S k (f,P) = k n I(f,P) + 0(k n - 2 ), (4.7) 
where the 0(k n ~ 2 ) error term stands for a distribution supported on P and bounded by a 
multiple of [|/[|can. For this, note first that if P is decomposed as a union of polytopes P\ and 
i"2 with disjoint interiors, then 

S k (f, P) = S k (f, Px) + S k (f, P 2 ) + 0{k n - 2 ) 

because the number of points of A k where there is a discrepancy is contained in the (n — 2)- 
skeleton of P\C\P2 and hence bounded by a multiple of k n ~ 2 . (The 0(k n ~ 2 ) error is also bounded 
by a multiple of supp(/).) In this situation we also have 

I(f,P)=L(f,P 1 ) + I(f,P 2 ). 

From these considerations, since we can decompose our polytope into integral simplices with 
disjoint interiors, it is enough to establish (4.7) for integral simplices. Although it is not hard 
to prove this by induction, we may simply invoke, for example, Theorem 1 of [KSW03] which, 
after rescaling, gives 

S k (f,Z) = k n I(f,Z) + k n - 2 E k (Z,f) (4.8) 
for any integer simplex S, where E k is a distribution on £ which is bounded by a multiple of 

ll/lbn. 

This is not quite the result we need, but if F is any facet of P, then we have 

S k (f,F)= Y f(a) + 0(k n - 2 ) (4.9) 
aeFnAfc 

(again because there are only 0(k n ~ 2 ) points in the (n — 2)-skeleton of F) and by what we've 
just proved, 

S k (f, F) = k n ^I(f, F) + 0(k n - 3 ). (4.10) 

Combining these observations with(4.7), we see that we can replace the sum over lattice points 
of the boundary by the corresponding integral, up to an allowable error term. 
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This completes our sketch proof. □ 

4.2. Divergence theorem. Apart from the Euler-Maclaurin formula, we also need a formula 
for the integral of the divergence of a vector field over the intersection of a hyperplane with P. 
In fact it is natural to consider a one-parameter family of parallel hyperplanes 

W(t) = {$(x)=t}, (4.11) 

where <J> is an affine-linear function on M n . In this situation we make the following definition: 

Definition 4.3. The number c £ M. is called a critical value of the one-parameter family P(~)W(t) 
if W(c) contains a vertex of P. If c is not a critical value of P n W(t), we call it a regular value 

of Pn w(t). 

We note that if to is not a critical value of PC\W(t), then for sufficiently small 5 > 0, PC\W(t) 
and PDW(s) are combinatorially identical and have the same conormals for s,t € (to - <5,to + <^)- 
It follows that if / is a smooth function on P, then 



/ fda t (4.12) 
Jpnw(t) 



t i-> 

ipnw(t) 

(which is continuous for all t) is smooth for t G (to — 6, to + 5). Here da t is the Leray form of 
PDW(t), i.e. 

dcr t d<£ = dx 

along W(t). 

Lemma 4.4. Let P be a convex polytope in W l and W{t) be as above, and suppose that to is a 
regular value of this one-parameter family. Let £ be a smooth vector field on P. Let P(t) be the 
part of P cut off by the half-space {&(x) ^ t}. 

Denote by da the Leray form of the codimension- 1 part of the boundary of P{t) and by dr 
the Leray form of the codimension- 2 part of the boundary. Then for all t in a sufficiently small 
neighbourhood of to, 

[ div(0 da = A / <£, d<&) da - [ (£, u)dr. (4.13) 

JW(t) at JW(t) JdW(t) 

Proof. Fix t and t + h near to so that the interval [t, t + h] contains no critical value (we assume 
h > here). Define the polytope C(h) to be the closure of P(t) \ P(t + h). The facets of C(h) 
are the two parallel facets P n W(t) and P n W(t + h) together with {C(h) n G}, where G is 
a facet of P. Denote by Z(h) the union of these 'side' facets of C(h). Applying the divergence 
theorem to C(h), we have 

/ div(0 d/x n = - / (£,d$)da t + / fodtydvt+h- f (4.14) 

We will now calculate the limit as h — > of this equation. 
By definition of Leray form, for any smooth function on P, 

f fdx= f +H \ [ fda s ) ds, (4.15) 
JC(h) Jt \JPr\W(s) J 

Thus if h is small, we have 

I fdx = h[ fda t + 0{h 2 ). (4.16) 
Jc(h) Jpnw(t) 

Similarly, for each facet G of P meeting W(t), we have 

f fda = h[ fdr t + 0(h 2 ). (4.17) 

JGc\C(h) JGnw{t) 



so that 

>Z(h) JiHJ'i-\\\\l)) 



[ fda = h[ fdr t + 0(h 2 ). (4.18) 

Jz(h) Jd(PnW(t)) 
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Thus the LHS of (4.14) is equal to 



h [ div(f) da + 0(h 2 ) 
Jw 

while the first two terms on the RHS combine to give 

h± f (t,v)da + 0(h 2 

M Jw(t) 

Finally the integral over Z(h) is 



h 



dW 



(e,^)dr + 0(/i 2 



(4.19) 



(4.20) 



(4.21) 



Combining these three equations, dividing by h, and taking h to thus gives (4.13) as required. 

□ 

With these preliminaries we can now establish Theorem [TTT] 

4.3. Completion of Proof of Theorem 1.1 We now complete the proof of Theorem |1.1| by 
deriving the distributional formula (1.14). Note first that that formula is written on X rather 
than downstairs on P. It is clear that ptk is T n -invariant, so it is enough to obtain ( |1.14| ) for 
functions / of the form //*(/), where / G C°°(P). Since the volume of each fibre of p is (2ir) n , 



H*(f)u n /n\ = f»» / fdx. 



x 



(4.22) 



Thus, identifying p t k and atk with their respective push-downs to P, we see that ( 1.14| ) is 
equivalent to the formula 



"(/r/ + M/p/ /dl + ( "'' / + ^ ,/> ) (/eC °° (P>) (4 



{Ptk J) = fe' 

where Rk is an appropriate error term and 



23) 



(at,f)= [ f^-~f fm 2 9 da. 
JN(t) 1 dt JN(t) 



(4.24) 



In the remainder of this section we shall always think of ptk an d atk as distributions on P rather 
than on X. 

With these preliminaries understood, we just combine the Euler-Maclaurin formula (4.6) 
with the distributional expansion (3.24), getting 

i-n-i^ - - r, 1 



(ptkj) = k r 



f + -k 
J 2 



dPt 



fda+ I ( s -^d l d 1 (G^f)))+0(k n - 2 ). 



(4.25) 



Now, by the divergence theorem, 



dP t 



dj{G ij f)vi da 



N, 



djiG^Dvida 



dP+ 



dj(G ij f)ui da. (4.26) 



Consider the second term on the RHS, and more specifically a facet F of P. We may suppose 
that F is given by x\ = 0, so the conormal is e\ and da = dx2 ■ ■ ■ dx n . 



From (3.12), 



Gei = 2x iei + Xl ^S{D)e 1 = {2x 1 + 0(x\))e 1 + x 1 r ] {x), (4.27) 



say, where the vector r](x) is orthogonal to e%. Hence 

fGu = f(2 Xl + 0(a;?))ei + x^x) 

and so 

dw(fGu) = 2f + 0( Xl ). 



20 



FLORIAN T. POKORNY AND MICHAEL SINGER 



By a similar argument to that used to prove (3.3), this is true uniformly up to the boundary of 

m from this 

didj(G ij f) 



F and so the contribution from this facet to the integral over dP t + is 2 J p f da. Hence ( gjgj l 
simplifies to 



Pt 



dj{G*f)vi&r-2 



fda. 



dP? 



Now use Lemma 4.4 on the integral over N(t) to get 

d 



Nt 



dj{G^f)uida 



dt 



N(t) 



fG ij ViVj dot 



(4.28) 



(4.29) 



the contribution from the boundary of N(t) being zero. The reason for this is as follows. Each 
boundary facet of N(t) is of the form N(t) n F, where F is a facet of P. By the lemma, the 
integrand will be /G y i/jZ/j, where v 1 is the conormal to F . But we have seen that G l3 v'- = on 
F in the previous part. Hence the boundary contribution is zero. Combining these calculations, 
we arrive at 

Combining equations ( |4.25 ) and (4.30), we obtain 

fdx + \k n ~ l I fsdx 
JPt 



N(t) 



f\d$\ 2 q da t -2 



fda. 



(4.30) 



(Ptk:f) 



+ 



k n 
1 




/ fm 2 g da] +0(k n - 2 ). (4.31) 
J N (t) J 



Now the first two terms in the middle line combine to give J N ^ f da, completing the proof of 
(423). The bound 

(R k ,f)^C\\f\\ C n + 4 

follows from the bounds on the error terms in (3.24) and (4.6). 

4.4. Complete asymptotic expansion. We note that these methods yield a complete distri- 
butional asymptotic expansion for p^: 

Theorem 4.5. There exists a sequence £j of distributions on P such that for each N > 0, we 
have 

(Ptk,f) = k n {bO.^ + ^W/)}, ( 4 -32) 

where i?j\r+l,fc is a distribution on P satisfying 

(R N+ i,k,f)<C\\f\\ cNl 

for some N' depending upon N . 

Proof. We sketch the proof as we will not use the result in the rest of the paper. Fix a test- 
gave an asymptotic expansion of (|e 0i fc| 2 ,/) 



3.5 



function / and an integer N > 0. Proposition 
to order N, and we took care to note that all the coefficients as well as the error term depend 
smoothly upon a. On the other hand, from the results of for example |KSW03] or |GS07j . for 
any smooth function u(a) on Pt, we have an asymptotic expansion of the lattice sum 

N 

u(a) = k n Y,Ij{v)k~ j + k- N - l {S N+1<k ,u), (4.33) 
p t nA* k j=o 

where the Ij and <Sjv+i,fe are certain distributions on P with the error term SN+i,k satisfying 

(S N +i,k,u) < C\\f\\ cN > 

for some integer N' depending upon N. Because all the coefficients in our expansion of 
(|e Q ,fc| 2 , /) are smooth in a , we may substitute (3.40) into (4.33), getting an asymptotic expan- 
sion of the form (4.32) □ 
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4.5. Combinatorial interpretation. Finally, we give a combinatorial interpretation of the 
quantity of the sub- leading term if / = 1: 

Proposition 4.6. If t > is a regular value of the family Pt, then we have 

Vol(5P+) = jf sAx -\j t J N fml^t- (4.34) 

Proof. If we plug / = 1 into the above formula, we know that we get the leading coefficients 
Aq{€) and A\(t) as the leading coefficients. On the other hand, the dimension of the space of 
sections is the number of lattice points in Pt and this is approximated by the Euler-Maclaurin 
formula with / = 1. Comparing coefficients now gives the result, at least if t is rational. Since 
both sides are continuous in t, the result is true for all regular values t. □ 



5. More general partial density functions and toric K-stability 



In this section we want to generalize Theorem 1.1 to more general subspaces of defined by 



a general rational convex polytope Pt of P, obtaining in particular a distributional asymptotic 
expansion for the partial density function p t k defined in this situation. 

We shall then use this distributional expansion to prove that toric cscK implies toric K-stable 



in a sense to be explained below. These results appear as Theorems |5.2| and 5.6 

We begin with a careful discussion of what we shall call 'polytopes with moving facets'. We 
have already seen the simplest example where a 1-parameter family of polytopes Pt is defined 
by intersecting a given polytope P with a variable half-space {^(x) ^ t}. This is a polytope 
with a single moving facet. We must generalize this to allow for an intersection of P with an 
arbitrary finite collection of half-spaces {<& a (x) ^ t} where each of the $ a is an affine function 
of x. 



5.1. Polytopes with moving facets. Let P C W 1 be a convex integral polytope. Suppose 
given a finite collection <i> a (x), (a S A) of affine-linear functions, with rational coefficients. For 
each t G 1, define 

P t = Pnf]{$ a (i)^t}. (5.1) 

a 

We assume that Pq = P. For any given t, define At C A to be the subset of 'effective constraints', 
i.e. a £ A t if Pt H {& a (x) = t} is a facet of P. With the assumption Pq = P in force, it follows 
that A = 0. 

Because the $ a are rational, given a positive rational number t, there is a positive integer N 
such that for all integers k divisible by N, Pt is a lattice polytope for the rescaled lattice A* k 
(equivalently kPt is a lattice polytope for A* for such k). 

Definition 5.1. For rational t > and integers k > as in the previous paragraph, define 
Vtk C Vfc to be the span of the sections of L k corresponding to points in P t n A£. Similarly, 
given a choice of toric metric h on L, ptk is defined to be the partial density function for the 
subspace Vtk- 

In addition to the combinatorial data P and Pt, we now choose a toric metric j on I. In 
order to state our generalization of Theorem |l.l[ we need the following notation and definitions. 

• Write dPt = N(t) U dP^~ , where N(t) - the 'new part' of the boundary - is the union 
of those facets defined by <& a (x) = t for a £ At- The Leray form is denoted, as usual by 
da. 

• A positive measure dp supported on the (n — 2)-skeleton of N(t) is defined as follows: 
for any pair of facets N a (t) = {& a (x) = t} and Nt,(t) = {^(z) = t}, define 

dp ab = |d$ a - d$ fe |2 dr ab (5.2) 

on N a (t)r\Nb(t), and define dp to be equal to dp ab on the relative interior of A^ a (t)nA r ; ) (t). 

• The notion of a 'regular value' and 'critical value' of the family Pt defined below. 
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For the last of these, note that the set At is locally constant in t in general, but will jump at 
a finite number of values of t. These will be called the critical values of the family and will be 
denoted cf 

= c < ci < • • • < c m . (5.3) 

By definition Pt = for t < and t > c m . A value of t not equal to one of the Cj will be called 
'regular'. The significance of this notion is that if t, s & (cj-\,Cj) for some j, then the facets of 
Pt and P s have the same conormals and are combinatorially identical. In particular, the Leray 
forms depend smoothly upon t for t in any one of these intervals. 

Given these preliminaries, we can state our generalization of Theorem 1.1 as follows. 

Theorem 5.2. Let P and Pt, and ptk be defined as above. Let 

C t = P\Pt, 

so that P is decomposed into mutually disjoint subsets Ct, Nt and Dt = Pt\Nt. 

Then the partial density function pth associated to P(t) has the following properties: 

Ptk(x) = 0(fc-°°) ifxeCt 

and 

Ptk(x) = Pk(x) + Oik- 00 ) ifxe D t . 
Moreover, the O's are uniform if x moves in a compact subset respectively ofUt or Dt. 
Let f 6 C°°(P). Then, provided is t is a regular value of the family Pt, 



(Ptkj) = k n 



f + 



1 

2k 



D, 



sf+(a t ,f) +0 



where 



(atk,f) 



I 

JN(t) 



fda 



1 d 

2di 



JV(t) 



r 

JN{t) 



/dp 



(5.4) 

(5.5) 
(5.6) 

(5.7) 
(5.8) 



and where 0(l/k 2 ) denotes a distribution such that (Rk,f) ^ Ck 



Qn+4, . 



Proof. Equations (5.5) and (5.6) are established following exactly the same argument as for 



their counterparts in Proposition 4.1 



Moreover, the strategy for obtaining (5.7) is exactly the same as for (4.23): the only difference 
is that the calculation of 



Pi 



did 3 (G ij f)dx 



is more complicated. 



Indeed, the first step is the same and the analogue of (4.28) here is 

d i d j (G il f)dx = 



p(t) 



N(t) 



di{fG^d^)da-2 fda 



(5.9) 



(5.10) 



On the right-hand side, we have used an obvious shorthand: the first term should more properly 
be written as 



ditf&djQjdcTa. 



Consider a typical term 



N a (t) 



(5.11) 



(5.12) 



in this sum. We want to use Lemma 4.4 to simplify this integral. For this, let s > t and consider 
P(t) \ P[s), which we think of as a neighbourhood of N(t). Decompose this set as a union of 
polytopes C a (t, s); C a (t, s) is defined to be the convex hull of N a (t) and N a (s). Because t is a 
regular point of the family Pt, N a (t) and N a (s) are combinatorially identical and their boundary 
facets have the same conormals for all s sufficiently close to t. See Figure 1 for an illustration 
of this construction: N(t) is BCDEF, N(s) is HJILM; and the C a (t,s) here are BCJHB, 
CDIJC, DELID and EFMLE. 
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Figure 1 



We claim that we can apply Lemma 4.4 with P replaced by C a (t, s) and W(t) replaced by 
N a (t). For this to be the case, we need to know that the 'side faces' of C a (t, s) do not vary as t is 
varied. Now a typical side facet is either the intersection of C a (t, s) with an old boundary facet 
F of dP + or else the intersection Z a b(t, s) with Cb(t, s) for some b ^ a. It is clear that F does 
not move with t. As for Z a b(t, s), any point x on it satisfies & a (x) = &b( x ) = t' i f° r t ^ t' ^ s. 
In particular, the hyperplane containing the facet Z a b(t,s) is given by <&b(%) — ^a(x) = and 
its inward conormal is d&b — d<3? a - So such a facet also does not move with t. 

Thus we can apply Lemma |4.4| to obtain 

/ div(£)da =^ / (^d<S> a }da a - [ {Mdr (5.13) 

JN a (t) at JN a (t) JdN a (t)ndP+ 

J2 I <£,d($ 6 -d$ a )dT a6 (5.14) 

b JN a (t)nN b (t) 



for any smooth vector field ^ on P. 

In the particular case that £ = fGdQ a , the integral over dN a (t) ndP^~ vanishes for the usual 
reason that Gv = on any facet of dP + with conormal v (cf. ( |3.3[ )). Thus we obtain 

/ div(/Gd$ a ) da a = A / \d<S> a \ 2 g da a - V / (/(d$ a , d($ 6 - $ a )) s dr a6 . 

•/iVa(t) dt ^AT a (i) ^ JN a (t)nN b (t) 

(5.15) 

Summing over a, we get the formula 

f div(/Gd$) da = A / Id^da + V /" /|d($ 6 -$ a )|2 d r a& . (5.16) 

The last term on the RHS here is /^y^ / dp by definition, so by combining this equation with 



(5.10), we obtain 

/ d i d j (fG«)dx = -~ [ fm 2 g da-2 [ /da-/ f dp. (5.17) 

JPt M JdPt JdPi JN(t) 



Substitution of this into (4.25) gives 



(Ptk,f) = k n [ fdx + \k n - 1 [ s dx + h n - l {a u f) + k n - 2 E k [fl (5.18) 
JPt 1 JPt 1 

(at,f)= [ fda-~ [ f\d$\ 2 g dv- [ fdp (5.19) 

JN(t) 1 dt JN(t) JN(t) 

as required. □ 



where 
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Setting / = 1, we obtain the analogue of Proposition 4.6 in this case: 



Proposition 5.3. Let g be any (smooth) toric Kdhler metric on X and let s = s(g) be the 
scalar curvature. Denote by Apt the above measure. If t is not a critical value of the family 
{Pt}, 

Vo\(dP+)= [ s-^-f \u t \ 2 da t - ( d Pt , (5.20) 
Jp t dt J Nt J Nt 

where s is the scalar curvature and v% is the conormal to Nt- 

5.2. Test configurations and K-stability. In |Don02j . Donaldson proposed a definition of te- 
stability for polarized varieties (X, L) and, in the toric case, related K-stability to boundedness 
properties of the Mabuchi energy. We shall not reproduce the exact definition here. The 
rough idea is to consider 'degenerations' of (X, L) to a (possibly very singular) polarized variety 
(Xo,Lq) with a C x -action. In this situation one can define the Donaldson-Futaki invariant F± 
of (Xq, Lq); then (X, L) is K-stable if F\ < for all possible degenerations. 

In the toric case, there is a subclass of toric degenerations which can be defined combinatori- 
ally as follows. Let the polarized toric variety (X, L) correspond to the convex integral polytope 
P. Now, given the data of the previous section, define $> a (x,t) = Q a (x) — t and consider the 
polytope r C (the last variable being t), 



r = P x [0,oo) np|{$(s,t) ^ 0}. (5.21) 



It is convenient to augment the defining equations for T by explicitly including &o(x,t) = t 
which defines the base t = of V. 

We refer to any such T as (the polytope corresponding to) a toric test configuration for (X, L). 
We note that the 'roof of T (see Figure 2(a)) is a union of n-dimensional convex polytopes N a . 
Then Xq in this case is obtained by gluing together the toric varieties corresponding to the iV a 
to obtain a singular variety. 

By definition, a product configuration arises when r is P x [0, oo) cut off (possibly obliquely) 
by a single affine function <&(x,t) ^ — see Figure 2(b). A product configuration is called trivial 
if the roof is horizontal, i.e. given by P x {c} for some c > 0. 




(a) A polytope defining a toric (b) A polytope corresponding to a 

test configuration and a non- product configuration, 

critical level set of Pt- 

Figure 2 



In [Don02j, the following combinatorial description of the Donaldson-Futaki invariant was 
given: 
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Proposition 5.4. Let (X, L) be a toric variety with moment polytope P and let V be a polytope 
defining a toric test configuration for (X,L). Then the Donaldson-Futaki invariant ofT is the 
coefficient of k~ l in the asymptotic expansion ofwk/kdk, where 

w k = N(kT)- N(kP), d k = N{kP) (5.22) 

and N(A) denotes the number of lattice points in the integral polytope A. 

The significance of this definition in relation to K-stability is as follows: 

Definition 5.5. Let (X,L) be a toric variety with moment polytope P. We say that (X,L) 
is K-polystable with respect to toric test configurations if for every toric test configuration 
corresponding to an (n + l)-dimensional polytope T as above, the Donaldson-Futaki invariant 
Fi is ^ 0, with equality if and only if T corresponds to a product test configuration. 

The reader is referred to [Don02, §4.2] for the details. 

Our next theorem gives a formula for the Donaldson-Futaki invariant, given a choice of toric 



Kahler metric g, which is very analogous to the formula for the slope given in Theorem 1.3 



Clearly T is closely related to the family Pf of polytopes (5.1). More precisely, let it denote 



the restriction to T of the projection (x,t) i— > t. Then 7r 1 (t) = Pt, where Pt is as in (5.1). By 



analogy with £5.1 let us introduce the following notation: 

• Write dT = N U dT + , where N is the 'roof of T, that is the union of the facets defined 
by the hyperplanes & a ( x , t) = and T + is the 'vertical part' of dT — the union of facets 
contained in sets of the form FxK, where F is a facet of P. 

Define a positive measure dp with support on the (n — l)-skeleton of N as follows: if 



N ab = T n {$ a = 0} n {$ b = 0}, define 

dpab = |d$ a - d& b \ 2 g df ab on the interior of N ab . (5.23) 

This makes sense because d$ a — d<& b = d& a — d& b is a 1-form on P; its length can thus 
be measured with the toric metric g. Then define dp on N to be equal to dp ab on the 
relative interior of N ab for all a ^ b. 
• We say that c is a critical value of n if 7r _1 (c) contains a vertex of T. It is easily seen 
that c is a critical value of ir if and only if it is one of the c,- of (5.3). 

Now define 

1 



*-wwi t « (5 ' 24) 

the integral being over the roof N of T. Since dp is a non- negative measure, 

A(r) ^ for any T (5.25) 

and 

A(r) = if and only if the roof of T has no codimension-2 faces. (5.26) 

In other words, A(r) ^ with equality if and only if T corresponds to a product test configu- 
ration. 

Theorem 5.6. Let (X,L) be a smooth polarized toric variety with moment polytope P C M. n . 
Let r C W l+1 be a polytope defining a toric test configuration for (X, L) . Then, for any choice 
of toric Kahler metric g in the Kahler class c\(L) on X, the Donaldson-Futaki invariant ofT 
is given by 

Fx = (Av r (pr^G?)) " Av P (s(g)) - A(T)) , (5.27) 

where s(g) is the scalar curvature of g, Av^(/) denotes average value of the function f over the 
set A and pr 1 denotes the vertical projection T — > P. 



The following is a simple consequence of (5.27): 
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Corollary 5.7 (cf. [ZZ08J). Suppose that X admits a toric cscK metric in the Kahler class 
c\ (L) . Then the Donaldson-Futaki invariant of any toric test configuration with polytope V is 
^ 0, with equality if and only if V is a product configuration. In other words the existence of a 
toric cscK metric implies that (X, L) is K-polystable with respect to toric test configurations. 

Remark 6. We follow the sign convention of |Don 02j rather than [RTQ6] here, so for us negative 
Donaldson-Futaki invariant corresponds to stability. 

Proof. (Of the corollary.) If the scalar curvature s is constant, then it is equal to its average 
over P and also to the average of pr*(s) over V. Thus these averages cancel from (5.27), leaving 

voi(r) 

F --2YoUP) A{T) - 



The result now follows from (5.25) and (5.26). □ 
5.3. Computation of the Donaldson Futaki invariant. We use the following observation: 
Lemma 5.8 ([Don02]). The large-k expansion of Wk is given by 

Wk = Vol(r)/t n+1 + - Xol(dT + )k n + Oik" 1 - 1 ). (5.28) 

Given this result, the main problem is to understand Vol(<9r + ) in terms related to the metric. 
Since the intersection of dT + with a horizontal slice P% is the 'old part' dP+ of the boundary of 
Pt, we have 

Vol(ar+)= / Vol{dP t + )dt. (5.29) 
Jo 



where c m is the largest critical value of 7r as in (5.3). Combining this with Proposition |5.3[ we 
shall obtain the formula: 

Proposition 5.9. With the above definitions and notation, we have 

Vol(dr+) = J prj(a) - Vol(T)A(r), (5.30) 
where pr x is the restriction to T of the projection (x,t) h-> x. 



Proof. We integrate (5.20) over each interval (cj_i, a,) and sum over j, getting 

Vol(dr+) = f wl(s)+ [ \v\ 2 da- V \\ \v\ 2 da- [ \u\ 2 da 




N(tJ) JN(t+) 



>N(t-) 

Now the integral over iVo is zero because — by definition — Nq is empty. Thus we have 



(5.31) 



Vol(ar+)- / pr* 1 (s) = - V / \u\ 2 da- / H 2 da) - / \v\ 2 da - / dp t dt. 
Jr j=1 \JN(tj) JN(t+) J JN(t-) Jo 

(5.32) 

The next lemma matches up the terms on the RHS of this equation with the terms in the 
sum defining A(r), cf. (5.24). Recall that the roof N of T is a union of facets -/V a and its 
(n — l)-skeleton consists of intersections of the form iV a & = N a n Nf,. We call N a f, horizontal if 
it is contained in a horizontal slice Pt and non-horizontal otherwise. The point is that if N a b is 
horizontal, then it has to appear as a facet contained in N(t) C dP(t) and moreover t has to be 
a critical value of it, since N a b C dPt certainly implies that Pt contains a vertex of V. On the 
other hand, if N a b is not horizontal, then it meets Pt non-trivially for t in some interval /, and 
for each t in the interior of /, Pt n N a b = N a b(t) is part of the (n — 2)-skeleton of N(t). These 
rather simple observations may nonetheless help with the following: 
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(a) A critical slice P c in the case (b) A critical slice P c in the case 

where dat and vt vary continu- where dat and vt do not vary 

ously through t = c. continuously through t = c. 

Figure 3 



Lemma 5.10. We have 




N(cJ) 



\d$\ 2 g da 



\d%da + 



N(c-) 



|d^>|2 dcr 



E 



N ab horizontal" ' Nab 

dpab- 



dp ab . (5.33) 



and 



d Pt dt= 

N„ h not horizontal 



(5.34) 



Proof. (See Fig. 3). Pick Cj < c m , consider any facet F, say, of P c and in particular the 
contribution F makes to the sum on the LHS of (5.33). Suppose first that F is the intersection 
of P Cj with a single facet N a of T. Then Ft = N a n Pt for t near Cj. Hence the Leray form and 
conormal of Ft vary continously for t near Cj. So such facets contribute nothing to the sum on 
the LHS of ( pT33ft . 

The other case to consider is that F = N a n N b . Note that this is necessarily a horizontal 
face of N. Suppose that a and b are ordered so that F t = Ptd N a for t < Cj but Ft = Pt D iV& 
for t > Cj (it always being assumed that \t — Cj \ is small) . Then the contribution to the sum on 
the LHS of (15331) is 



/ (\d$ a \ 2 g da a -\d<t> b \ 2 g da b ). 

J F 



We claim that this is equal to 



N a nN b 



|*a - $b\g df ab 



dpab- 



Because N a and A^;, meet in a horizontal plane, we can choose affine coordinates on 
$ a (x, t) = Axi - t, 5>b(x, t) = Bx\ - t, with A > B > 0. 



(5.35) 

(5.36) 

so that 
(5.37) 



(Thus r is given locally by the intersection of the half-spaces Ax\ ^ t and Bx\ ^ t.) Then 



d$ a = -Adxi,d$ b 



1 1 

-Bdxi, da a = -rdx 2 ■ . . dx n , da b = —dx 2 ■ ■ . dx n , 

A ID 



and so the integrand in (5.35) is 



(A - B)\dxi\ 2 dx2 ...dx 



(5.38) 
(5.39) 
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On the other hand, 



|d$ a -d$ 6 | 2 = (A-5) 2 |dxi|2 df, 



1 



ab 



-dx 2 ■ ■ . dx n , 



(5.40) 



A-B 

from which the equality of (5.35) and ( |5.36| ) follows. 

This proves (5.33) except that we have ignored the last term on the LHS and the codimension- 
2 faces of T contained in t = t m . The required equality can be obtained by straightforward 
modification of the above discussion and further details are omitted. □ 



The proof of Proposition 5.9 follows by combining this Lemma with the definition of A(T) . □ 



We can now easily complete the proof of Theorem 5.6: just substitute (5.30) into (5.28) and 



use 



(5.41) 



kd k = Yol(P)k n+1 (l + ^ Avp^AT 1 + 0(AT 2 )^ 
to calculate 

W t = W) + 2W) iAM ^ * W) - AT ' W - A(F)) + ° {k ' 2) - 
The result now follows from the definition of F\ as the coefficient of k~ 1 in Wk/kdk- 
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